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ABSTRACT 


In recent years the subject of suppressing aerodynamic 
noise by shielding has actively been pursued both analytically 
and experimentally. Although there is a reasonably good 
agreement between the experimental results obtained under 
laboratory conditions and analytical results » the comparison 
is far from being satisfactory under flight conditions where 
several physical factors come into play in this problem. In 
order to make a reasonable study of the problem, the whole 
range of the frequency spectrum has to be studied. 

In the present study the frequency spectrum has 
conveniently been divided into two regimes j the low frequency 
and high frequency regimes. Two separate methods have been 
developed for application to each regime. 

For the long wave length propagation, the acoustic field 
due to a point source near a solid obstacle may be treated 
in terms of an inner region where the fluid motion is 
essentially incompressible, and an outer region which is a 


viii 



linear acoustic field generated by the hydrodynamic 
disturbances in the inner region. This method has been 
applied to a case of a finite slotted plate modelled to 
represent a wing extended flap for both stationary and 
uniformly moving medium. 

In the case of short wave length propagation, a very 
effective approach utilizing a combination of the method of 
ray acoustics, the Kirchhoff’s integral formulation and the 
stationary phase . approximation has been developed. The 
examples studied using this method include many limiting 
cases. The solutions of these limiting cases agree with the 
known solutions. This method, too, has been applied to a 
new problem of physical interest. The problem consists of a 
semi-infinite plate in a uniform flow velocity with a point 
source above the plate and embedded in a different flow 
velocity to simulate an engine exhaust jet stream surrounding 


the source. 
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Chapter I 
Introduction; 

...."Nature operates by the simplest and 
most expeditious means".... 

-Fermat- 

The introduction of more powerful jet engines in civil 
air transportation in the 1950s drastically revolutionalized 
the aircraft operational speeds and altitudes. The side 
effect of this advanced technological advance was the build 
up of noise on approach to land and take off paths on 
airports. 

Although the later generation of engines (turbofans) 
are relatively quieter, the increased sizes and power 
delivery of such engines still caused unacceptably high 
sound pressure levels (generally of the order of 110 dB at 
30.5 meters). With a desire to operate aircraft services 
close to the city centres and to use shorter runways, the 
Short Take Off and Landing (STOL) concept using enhanced 


lift produced from the interaction of the engine jet with 
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the wing and flaps has evolved. This arrangement is used 
now to shield engine noise from an observer. The concept 
has attracted a wide range of research activities both 
analytical and experimental. 

The numerous studies on diffraction and noise reduction 
by barriers can be divided into three groups roughly; 
theoretical studies, experimental studies using scaled 
models and full scale experiments under normal enviromental 
conditions . 

Comparison of experimental (under laboratory 
conditions) and theoretical results shows some disagreement 
for the lower frequency although there is a good agreement 
for the high frequency end of the spectrum. 

A full scale experimental study was extensively carried 
by Jeffrey and Holbeche (19) using a Delta wing aircraft. In 
the long v;ave length region, particularly, the results 
showed a marked departure from those obtained by scaled 
model experiments or the theoretical ones. Similar 


observation is seen in experiments by Conticelli, De Blasi 
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and O'Keefe (7). Hoch (16) used the Bertin Aerotrain to 
study the forward flight effects on aircraft engine noise. 

The results are in agreement with other studies although, 
again, there is a disagreement for the low frequency 
spectrum. Fink (15) reports results that show similar 
phenomena. 

There are many physical factors that exist in the real 
situation but are generally not included in the analytical 
studies. In the following discussion, a few of the more 
important ones are mentioned. 

1. Sound Source ; In the analytical studies, the 
source is, in general, simulated mathematically by acoustic 
poles which are also usually assumed concentrated at a point. 
In the real situation, the sound is generated aerodynamically 
over a certain region by the mixing of the high speed 
turbulent jet with the free stream. Internally, the 
rotating compressor and the turbine blades are also sources 
of noise v/hich transmits outside the engine. The 
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roechanisiTis of generation and transmission in the source 
region are extremely complex and by and large of an unknown 
nature. The strength, frequency content and directivity of 
the generated sound are not possible to be determined 
analytically to any degree of accuracy by the present 
knowledge. 

2. Transmission Path: Before reaching the atmosphere, 

the sound generated at the source region must transmit 
through the jet exhaust. Due to the temperature, density 
and velocity non-uniformity and the temporal variations of 
these quantities caused by turbulence, the transmission 
characteristics of noise in this "near field" region is not 
easy to determine theoretically. 

3 . Shrouding Effect : At the interface between the jet 

and the free stream, there are strong temperature and 
velocity gradients resulting in reflection and refraction 
of the sound v;aves and causing the intensity, directivity 
and frequency content to change in the far field from those 
pfg^icted by any simple theory assuming a uniform atmospheric 



condition. 
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4 . Reflection and other Effects due to the Geometry 
of the Wing and the Airframe ; For the engine 

under the wing configuration, for example, sound may be 
strongly reflected by the under surfaces of the wing and the 
flaps. This causes a large increase of the sound intensity 
received on the ground. For the engine over the wing 
configuration, on the other hand, noise is greatly reduced. 
During landing, the gap between the main airfoil and the 
extended flap may cause "leakage" of sound from the upper 
side of the wing to the under side of the wing and thereby 
reducing the effectiveness of wing shielding. 

5. Frequency on Diffraction ; The aircraft noise 
spectrum extends from very high frequency sounds to very low 
frequency ones. The low frequency sounds are not very 
effectively shielded by the wing because of the strong 
diffraction of waves of long wave length. 

6. Forward Flight Effect ; To compare laboratory and 
in-flight data, it is important to know how the motion of 
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the aircraft affects the predicted or measured results. 

7. Effect of Shed Vortex: Earlier investigations on 

wave-barrier interaction did not make any note of the sound 
that can be generated by shed vortex for finite or . 
semi-infinite surfaces. This issue has occupied several 
reseachers in the last twenty years: Orzag et al (35) ; 

Jones (20); Davis (11); Jeffery et al (18); and 

Broadbent (2) , among others. 

These studies have indicated that the trailing edge 
shed vortex has three roles; 

(i) It acts as a shield to sound from one region of 
such a vortex to the shadow. 

(ii) Such vortex is itself a source of sound. 

(iii) The vortex does refract sound if the source of 
sound is close to the vortex, but it has less 
refractive effect if the source is far enough, 
Cooke (8) . 

8. Effect of Viscous Boundary Layer; Powell (36) 
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studied the problem of sound caused by boundary layer on the 
solid surface. This is what many experimental investigators 
have referred to as scrubbing noise. In this case the 
dipole source noise as well as the quadrupole sources are 
present. 

Many of the above phenomena are studied in classical 
theories in the propagation of light waves. The study of 
light refraction, ref laction and diffraction has a long 
history. In the field of optics, the basic principles 
enunciated by Fermat, Huygens and Fresnel have been 
formulated into a general integral theorem by Kirchhoff. 
These basic principles, borrowed from the field of optics 
and later generalized in the general electromagnetic wave 
propagation after the Maxwell's development of the field 
equations , have been found applicable in the study of short 
wave length sound propagation. A combination of these 
principles has been employed in formulating a method 
applicable for short wave length acoustic propagation in the 


present study. 
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In recent times, with the rapid expansion of thO air 
transportation, the subject of acoustic wave •> interaction 
with a barrier has attracted a considerable amount of 
research. Various theories have been explored in the 
ensuing search for various solutions to the problem. 

The various published theories differ mainly in the 
way in which the fluctuating flow is assumed to interact 
with a barrier (with an edge of a plate in particular) to 
produce sound and fall roughly into the following four 
categories: 

(i) General theories based on aerodynamic noise based 
on the Lighthill (1952) acoustic analogy. These 
have been developed by Curie (10); Ffowcs-Williams 
and Hawkings (14); Howe (17); Jones (20); 

Powell (36); Ribner (39) among others. 

(ii) Theories based on the solution of special 
problems have been tackled by: Broadbent (2,4) ; 


Cooke (9); MacDonald (31); Orzag and Crow (35). 
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(iii) Theories based on ad hoc models involve postulation 
of source distribution whose strengths and 
multipole types are generally determined 
empirically. This group would include: Clapper et 
al (6); Lan et al (22); Larson et al (23) and 
Maekawa (32) among others. 

(iv) The experimental work that generally may include 
analytical background material may be found in the 
work of Conticelli et al (7); Fink (15); Hoch (16); 
Jeffrey et al (18); Lush (30); Reshotko et al (38) 
and Strout and Atencio (40) . 

The above list is no way exhaustive but it is 
representative of the work that has been performed in 
different approaches. In all cases, however, the main 
concentration has been the study of short wave phenomena. 

The present study, presented in chapter II takes the long 
wave length into account to bridge the apparent gap in the 
literature. 

In the chapters to follow, the wing shielding problem 
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will be treated in two aspects: One aspect involves sound 

of long wave length and the second aspects deals with sound 
of short wave length. 

The main purpose of the analysis is to develop 
mathematical methods and to find physical solutions for the 
wave-barrier interaction problems typified by Sommerfeld's 
classical treatment of diffraction by a semi-infinite plate. 
The particular interest here is in examining the effectiveness 
of using the wing as a barrier to shield noise from the 
engine installed above the wing under flight conditions. 

In chapter II, the shielding of a concentrated 
quadrupole sound source by a semi-infinite plate with a 
narrow slot is considered. This problem simulates the 
shielding by a wing with the flap in the extended position. 
The acoustic field as affected by the slot is determined by 
the method of matched asymptotic expansions, under the 
assumption of long v/ave length approximation. 

In the outer region, far av;ay from the slot, the 
presence of the slot can be neglected for a first 
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approximation. In such a case the solution is easily 
obtained from the principle of ray acoustics. Such a 
solution, however, would not be satisfactory in the region 
close to the slot. 

The solution close to the slot is calculated by close 
examination of the inner flow. The long wave approximation 
together with the Prandtl-Glauert transformation enables one 
to derive the inner solution from the Laplace's equation. It 
is found that the inner solution contains a second order 
term which is a dipole singularity. In order to match the 
inner solution to the outer one, it is necessary to calculate 
the second order terms for the outer expansion. The dipole 
source term in the inner region becomes the generator of the 
secondary acoustic disturbances in the far region. It is 
observed, in conclusion, that the location of the engine 
noise source in relation to the slot is crucial in the 
shielding effect of the noise by the wing with a flap. 

A method, suitable for deriving solutions for cases of 
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short wave length, is developed in chapter III, The method 
combines the basic techniques of ray acoustics, the 
Kirchhoff's integral and the stationary phase approximation. 
The basic principle of the method is based on the problem of 
a semi-infinite plate with a sinusoidal monopole source 
located above it. In the limit of very large wave number 


, the approximate solution can be divided into three 
regions with sharp discontinuities. The solution is 
reasonably accurate and the order of the neglected term is 


0(-^) . This solution is not accurate at the boundaries 
separating the three regions where a transition layer with 

/ i X 

contributions to the solution to the order of are 


present , 


The method proposed for use is, therefore, to use the 


Kirchhoff's integral with a surface of integration lying 


outside the transition layers. One can then use the 
solution obtained by ray acoustics for Cp and 

Kirchhoff's integral to determine the function cj) at any point 


P, including that within the transition layer. The integral 
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generates the second order correction term which is of order 
. The method of stationary phase is used to 
determine the second order approximation directly from the 
Kirchhoff's integral. In a sense, the method is to treat 
the Kirchhoff's integral as an integral equation and then 
solve it by the iteration procedure, using the ray acoustics 
as a stating first order approximation. 

It is found that the stationary point "v/indow" is of 
the order of where ^ is the wave number of the wave. 

The acoustic wave received at a point will depend on the 


proximitly of the stationary point window to the trailing 
edge. It is then determined that the solution at a point 
P (observer position) is basically that due to the ray 
acoustics modified by the shielding factor F (^) ♦ The 
shielding factor is a Fresnel integral form and cuts down 
the amount of the intensity of the wave as the stationary 
point window moves close to the edge, so that in the limit, 
the edge blocks the stationary point window leading to the 
total zone of silence (shadow) . 
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The method is tested by solving two classical problems 
of wave diffraction. The first one is the interaction of 
the sound wave from a three dimensional point source with a 
semi-infinite plate first studied by MacDonald (31) and 
later extended by Cooke (9) under the condition of large k » 
which is the same as the condition of large distance used 
in the present study. The result obtained is in complete 
agreement with equation 13 on page 9 of Cooke (9). The 
second problem is the interaction of a plane wave with a 
semi-infinite plate. The far field solution obtained by the 
present method is in agreement with equation 2.86 on page 73 
in the book by Noble (34) and with the result obtained more 
recently by Candel (5) . 

The method developed in chapter III has been applied to 
solve a problem of practical interest in chapters IV. The 
problem involves two flows with a semi-infinite plate 
separating the two uniform flows of different velocities. 

A source is assumed imbedded in the flow over the upper 
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surface of the plate. This is intended to simulate the 
problem of an engine situated above the wing. The acoustic 
field in the region below the wing, caused by the 
transmission and diffracted fields is to be determined. The 
boundary conditions on the free interface regarding the 
continuity of pressure and particle displacements determine 
the reflection and transmission coefficients of the wave at 
the interface. Using the Kirchhoff's integral for the 
semi-circular boundary and the theory developed in chapter 
III, the solution is derived without much difficulty. It is 
found that, in addition to the shielding factor F (p) 
dependence, such a solution is dependent on the reflection 
and refraction of the incident ray at the interface. The 
Mach number has the effect of changing the boundaries and 
the thickness of the transition zone. 

The main conclusions and applications of the developed 
principles are outlined in chapter V. From the results 
obtained in the problem of a semi-infinite plate, results 
for a finite plate as a model for an aircraft wing are 
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derived and numerical results presented. 

The effect of the gaps, angular deflection and other 
geometrical variations of the wing surface, which represents 
the real physical situations has been studied and it is 
observed that the effects of gaps on the far field is 
equivalent to that produced by localized acoustic poles 
induced at the position of the slot. 

Another major interesting factor regarding the effects 
of source distribution which has been studied by Thiessen 
(41); Levine (27); and Embleton (12) is re-examined for the 
geometry of the semi-infinite plane. The appropriate power 
radiation in the different regions determined by ray 
acoustics is presented in integral form. 

Finally, the question of how the flow velocity affects 
the shielding of a semi-infinite plate in a subsonic flow is 
briefly studied. It is found that the fluid motion causes 
two important effects: 

(i) The boundary of the shadowed region depends only on 
the location of the source relative to the plate 
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trailing edge. 

(ii) The angular thickness of the transition zone is 
affected by the fluid motion. The quantitative 
values of these effects can be numerically 
evaluated as a function of the Mach number , as 


shown by a numerical example. 
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Chapter II 

THE SLOTTED WING PROBLEM; 

II-1 Background ? 

In recent years the concept of the Short Take Off and 
Landing (STOL) has evolved a principle of externally blown 
flap to generate the additional high lift necessary for such 
operations. In the approach or take off phase, the jet from 
the engine interacts with the extended flaps. The slots 
between the main wing and the flaps become sources of noise 
as shown in figure (2-A) . The engine may be situated above 
or below the wing. Figure (2— A) shov/s the case of an engine 
below the wing but the principle applies for the over-the- 
wing extenal blov/ing as well. 



Figure (2-A): Under-the-wing jet blowing. 
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Although a substantial amount of effort has been made 
experimentally to determine the noise emitted from the engine 
jet and the wing/flap interactions (see Lasagna et al (24) 
and Falarski et al (13)) where it has been observed that the 
deflection of the flap generates a greater amount of sound 
intensity, there has been a limited analytical investigation 
of the problem. Ting (42) and Leppington (25) have 
investigated curvature effects of barriers on wave 
diffraction but such work did not involve gaps (slot) 
effects which is the real situation in the take off and 
landing phase of aircraft operation. The mathematical 
model, and method used in the study of this problem is 
presented in section II-2 below. 

II-2 Mathematical Model ; 

In order to analyse the problem described in section. 
II-1ja model is constructed as in figure (2-B) next page to 
represent a slotted semi-infinite plate and a quadrupole 
source S. A uniform flovv v/ith Mach number M<1 is assumed. 
The source S is monochromatic with frequency f = ^ 
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where X is the wave length and is the speed of the wave. 
The dimensions h, b, and L are assumed much larger than the 
wave length of the propagating wave X • ^be plate is 
assumed to be of zero thickness and zero angle of attack. 



Figure (2-B): Slotted wing and noise source in uniform 
subsonic flow. 


A method of matched asymptotic expansions using k£=:£ as 
the small parameter will be applied in the solution of this 
problem. In the region far from the slot ( r»,£ ), the 
presence of the slot may be neglected in the first 
approximation. The solution for this region is referred to 
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as the outer solution. For convenience/ the outer solution 
is constructed by the ray acoustics on the account that by 
assumption ^ h/ b, and L. As a result/ the solution for 
the outer region has discontinuities at the boundary of the 
regions according to ray acoustics and it is necessary 
therefore to relax condition (ii) in section II-3 to allow 
such discontinuities. In the small neighbourhood surrounding 
the slot/ the time variation of the fluid motion is nearly 
simultaneous and therefore a guasi-steady approximation may 
be applied to the wave eguation. The resulting eguation is 
the steady compressible flow equation which can be reduced 
to the form of the Laplace equation by the Prandtl-Glauert 
transformation. In solving for the inner solution near the 
slot/ the exact boundary condition of the slot must be used, 
however . 

In the following analysis the total velocity potential 
will be written as the superposition of the uniform flow 
potential Ux and a perturbation potential f / so that: 

^ = Ux + ^ 


( 2 - 1 ) 
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where ^ is the potential due to the presence of the source 
and the slotted plate. The solution of cf) is sought in the 
analysis that follows. 

II-3 Field Equations ; 

The governing equation is the convected wave equation: 


%+b)5Cy-h) 


(2-2) 


where <p is the velocity potential defined in equation (2-1) 


with 


^ J j m • i-u 1 -j. j c 

— g-.r: Kr"? 5 ^ A ' 1 -VT 7 - IS the amplitude of 

the source strength in tensor notation. The boundary 


conditions for this problem are 
(i) 


O on y = 0; —03 < X and 

^ V2 < X ( I* + V2 ) . 

(ii) Cp f ) <^^are continuous everywhere except on 

the plate surface. 


(iii) The Sommerfeld radiation condition requiring 


outgoing v;ave at infinity must be satisfied. 


II-4 The first order outer solution. 


The first order outer solution will be derived by the 
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use of ray acoustics approximation. (This is not a 
necessary assumption, but it will serve to give an explicit 
solution to illustrate the procedure better) . The outer 
region is, by definition, one which is far from the slot. 

In this region, one may assume an outer expansion of the 
form: 

0 = ^ + ef + t + 

The first order solution for 0 corresponds to the case 
when 0; that is, cj) can be obtained from equation (2-2) 

with the approximation that the slot does not exist. 

The solution for ® for a semi-infinite plate and a 
quadrupole source above it is derived by the method of ray 
acoustics. The acoustic field due to the source, according 
to this method, can be divided into three regions: a region 

where both the incident and reflected waves are present; 
another region where only the incident wave prevails and 
finally one finds a region (below the plate) which is in a 
complete shadow^ figure (2-C) next page. 
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— cx> 



Figure (2-C) : Sound emission from a source in the presence- 
of a semi-infinite plate. 


The solution for a monopole source at position ( -b,h ) 
in uniform flow is known in the literature ( ref. 37 ) , and 
is given as: 


ikMX r 

* - [ «: 


(2-4) 


where r = 


IV2 


(X + b) +(l-M^)(y_h) ana 


being the Hankel functions of the first kind and zeroth 


order . 

The solution for a guadrupole can constructed using 
equation (2-4). For the present purpose, only solutions of 
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(p in regions I and III are needed. In region I above the 
plate, is given by: 






(2-5) 




where r = |^(^X + b) + C X ^ 

In equation (2-5) above, the term involving ^ t-M^ y 

gives the potential due to the quadrupole in the absence of 

the plate; while the second term involving W i - — ^ ) is 

o >1 —My 

introduced to satisfy the solid boundary condition on the 
plate surface. 

The appropriate sign in equation (2-5) is chosen in 
such a way that for quadrupoles whose directional 

. a"* <3"' , . . 

derivatives involving -- — r or -r — pterins, a plus sign is 
selected between the two terms; while a minus sign is used 


for quadrupoles whose directional derivatives involve terms 
of the form... ■ — . 

3x3y 

According to ray acoustics, however, there is no 


acoustic perturbation in region III under the plate; so that: 
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$ = o 


( 2 - 6 ) 


From the first order solutions equations (2-5) and (2-6) for 


the acoustic perturbation/ the velocity field can be 


determined. The perturbed velocities for the upper and 


lower surfaces of the plate are given by: 

^ d ^ 

over the plate u = I — — * V = O 

“ \dX J ^ ° 

^ /y=o 


below the plate 11=0 3 V r= O 

o o 


(2-7) 


The velocity is, therefore, tangential to the plate and 
discontinuous across the plate. Such a discontinuity is 
acceptable across the plate but is not acceptable across the 
slot. Near the slot, a region referred here as the inner 
region (see figure 2-D) , an inner solution must be 
constructed . 

The velocity component at the slot is denoted as: 



X =0 
y=ro 


( 2 - 8 ) 
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where X - 0 and y = O is the location of the centre of the 
slot. U is an implicit function of the location, strength 
and direction of the quadrupole source. The velocity field 
in the inner region should approach IJ on the upper side of 
the plate, but will go to zero in the lower part of the 
plate. 

1 1 -5 The first order inner solution. 

The inner region is a region of the dimension O(^) 
containing the slot figure (2-D) . Since the wave length is 
much larger than the width of the slot (A »£ ),'.the time 
variation is nearly simultaneous for all points inside this 
region; So that the differential equation to the first order 
approximation is a steady flow equation. 



I 

Figure (2-D): Configuration of the slot region. 
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More precisely, one may write the potential (j> in terms of 
the following inner expansion: 

(p ^ ^ • (2-9) 

Substituting such an expansion into the governing eguation 
(2-2 ) 7 and after the quasi-steady state approximation is 
applied/ the first order inner solution (h satisfies the 
following Prandtl-Glauert equation: 



For the potential (h , the exact boundary conditions at the 

*0 

slot will be used. However, since the slot is narrow, the 
plate may be treated as infinite in both directions. In 
addition, must match the outer solution at large 
distances; so that the boundary conditions may be stated as 
follows ; 


(i) 




ay 


= O on y = O , -00<X<-^/2 and ^/2<x<^co 


. t d(f> 

continuous in the flow 


region. 


(iii) (j) is to match the outer solution. 
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This last condition requires that the velocity field 
obtained from ^ must match with the outer velocity field, 
thus: 


ii =0 

ax 


as y — > + CD 

as y — »-cD 

The above boundary value problem may be solved by 


and ^ O 

ax 


introducing the complex potential ^ ^ given by: 

where ^ ) and ^ ) are the velocity potential and 

stream functions respectively. The solution ( Z^ ) 


adapted from ref. 33 is given by: 


Fo ( Z 


, ' =u/2 [ z, . ] 


( 2 - 11 ) 


where H, = x,+iYi s x,= x , Y,=yyi-M* , 


The funct 


ion ^ ^ is defined as 


where , 6^ , 


i 5ti!r 

e 2 = a + i b 

and 0 are as shown in the figure 


(2-E) 
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below. The functions <f>^ and can now be written as: 



where 9^ and 0_ are such that the point P should not move 
across the cut in figure (2-F) next page. 

In terms of cp , the analysis above gives; at large 


distance from the slot: 



Figure (2-E): Geometry of the slot for large distance' 
solution calculation. 



Figure (2-F) T he definitio n of the complex function 
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The first term of the solution given by equations (2-12) 
matches with the velocity of the first order outer flow 
field cf) ; while the second term has a form of a dipole and 
represents the effect of the slot. 

The inner solution has some rather interesting features. 
In. equation (2-11), the first term y(Z represents a 

uniform flow of velocity U/z along the direction of the 



Figure (2-G) : Uniform flow. 



I 


Figure (2-H) : Flow past a slot. 
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The streamlines are adapted from (Ref. 33) with a change 
in the branch cut. When the uniform flow represented in 
figure (2-G) is superposed to figure (2-H ) , the flow 
velocity in the upper plane approaches the uniform velocity 
Uat large distances; while in the lower region, the flows 
in the two figures are in opposite direction and tend to 
cancel out, so that only the flow velocities of lower order 
are present. The resulting flow pattern is shown in figure 
(2-1) below. Attention is called to the fact that the flow 
depicted here represents the perturbation potential <j) of 
equation (2-1 ) . 


tx 
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It is noted that at large distances, the flow field 
under the plate looks like that of a dipole with a source at 
( - 0 ) and a sink at ( ^/2, 0 ) figure (2-J) below. 
This flow is represented by the first term in the asymptotic 
expansion in equation (2-12) (b) . In the upper half plane, 
the streamlines bend downward tov/ard the plate on the left 
side of the slot and turn away from the plate on the right 
side of the slot. At large distances, the flow pattern is 
equivalent to the superposition of a uniform flow of 
velocity U and a dipole with a sink located on the left edge 
( - ^/2, 0 ) of the slot, figure (2-K) . 


Figure (2-J) ; Flow under the plate 
in the slot region. 




Figure (2-K) : 


Flow over the plate 
in the slot region 
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Thus the far field generated by the slot in a basic 

/X/ 

flow of velocity t/ fory>0; and zero velocity fory^<0 is 
of a rather peculiar nature, namely, the perturbation at 
large distances are the pertubation given by dipoles of 
different signs. Note also that at both edges of the gap 
( X|~ + '^/2, 0 ) the X -- velocity component has a well known 
square-root singularities which arise due to the assumption 
of the plate being of zero thickness. 

On the line y^s 0, downstream of the plate, the 
velocity fields generated by the tv;o dipoles of opposite 
signs are discontinuous. The velocity components given by 
the two dipoles are both zero at y^ = O and hence are 
continuous, but the U— components are non-zero and of 
opposite signs and hence discontinuous. Since — O is the 
boundary between the acoustic field generated by the two 
dipoles ( one on the upper surface and one on the lower 
surface of the plate ) this discountinuity is acceptable 
according to ray acoustics. 

There are three important observations to be made about 
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the inner solution. 

(i) The slot acts like a dipole. 

(ii) The axis of the dipole is along the surface of the 
plate (X-axis) . 

(iii) The signs of the dipole are different for the 
region y>0 and y<0. 

A simple explanation of these features is given as 
follows; 

(a) Due to the long wave length assumption ( ^/^«d ), 
the effect of the slot is mathematically equivalent to the 
scattering of a plane wave by the incompressible flow region 
surrounding a small particle (Rayleigh scattering of scalar 
waves). The inner solution surrounding the slot, therefore, 
acts like a dipole to the outer acoustic field. 

(b) By the assumption of a thin plate lying on the Xpaxis> 

the dipole must lie on the X-axis ( y s= O ) with an axis 

‘ I 

directly along the plate. A vertical dipole would generate 
a vertical velocity component at = O and this cannot 
satisfy the boundary conditions imposed on jr = O, 
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(c) The reason that the dipole has a different sign is 


mathematically due to the choice of the branch cut for the 


function^Z^ — £ A as shown on page 31 . 


In the branch cut as in figure ( 2 -L) below, the 


infinity point will not be contained in the cut. In such a 
case, the real part of Zj — £ /4 denoted by Q will 


approach the same value X| both for » 00 and y^- 


-00 • 



Figure ( 2 -L); Infinity point not in the cut. 


In the present problem, however, the matching condition 
imposed by the outer solution requires that C|-^ Xj as ^00 
and O — ►— Xj^as y^— ♦>— 00. in order to obtain this result it is 
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necessary to draw the branch cut from to -oo and 

from to + 00 . in this way, the infinity points of 

the two upper and lower planes, in figure (2-F) are on two 
different Riemann sheets in figure (2-L) so that the 
behaviour of the function and hence the dipoles have different 

sign for and y<0. Physically, the reason that the 

1 * 


matching conditions reguire the real part of to 

behave like X^for y^>0 and -Xj^ for ^<0 is that the outer 


solution for the upper and lower regions are obtained by ray 
acoustics, and there are strong discontinuities at the 
interfaces of the different regions. 

11-6 The Second Order Outer Solution ; 

The differential equation for the second order outer 
solution is: 






dy' 



= O 


(2-13) 


The appropriate boundary conditions are: 


(i) 


(p ) ^ u must be continuous except 

' dx 'dy 
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at the plate and at the boundaries according to ray 
acoustics. 


(ii) Radiation at infinity must be satisfied. 

(iii) (f) must match with the inner solution. 

This means that as X and y approach the position of the slot 
( X = y = O ) ; takes the asymptotic values below. One 
must note that the discontinuity of the solution on y = O 
for X>£/2 must be admitted just like the discontinuities on 



follows : 

Considering in region I above the plate^ one may write: 
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Where /\ is a constant to be determined. The above 
expression for 9 ^ represents a two-dimensional acoustic 
dipole in a uniform flow, and is obviously a solution of the 
governing equation (2-13). Such a solution also satisfies 
the radiation condition at infinity since it represents an 
outgoing wave. To determine the constant A by the marching 
condition one must expand the expression for small values 
of X and y to obtain: 


^ ax 


^ kgAi 
T ax 

X 




(2-16) 


K 


x'" + (i-M ) y 


Comparing expression (2-16) above with the inner solution 
as presented in equation ( 2 - 12 ) one obtains: 

zktA^ _ 

T 16 


SO that 


A - 


irU 

3Z k 


(2-17) 
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Hence 




(2-18) (a) 


Similarly, one finds that in region III below the plate: 



irfUi ± 

3 ^ dy^ 



(2-18). (b) 


I 1-7 Discus s ion and Conclusions : 

(1) From the equations (2-18), one observes that the 

effect of the slot is equivalent to a dipole situated at the 
position of the slot. The strength of the dipole is 
directly proportional to U which is the velocity at the 
surface location of the slot due to the quadrupole. U is a 
scalar function of the strength, directivity and location of 
the quadrupole T.. , and may be determined by the use of 
equations (2-5) and (2-8) , when the value of is 

explicitly given. Note that the dipole is of opposite sign 
on the two sides of the slot figure (2-1) . 

(2) The sound field pattern generated by the slot in the 
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uniform moving medium is sketched in figure (2-M) below and 
the relative amplitude pattern for the velocity potential is 
in figure (2-N) below, 

(3) The strength of the induced dipole is proportional 
to { Ji^ ), so that if the slot width Z is reduced by half 
(say), the strength of the dipole is reduced four times. 



Figure (2-N); Relative amplitude of the velocity 
potential |9l* 
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(4) The asymptotic behaviour of the Hankel function of 


argument W is: 



as Of —►00 


Hence the strength of the induced dipole is proportional 

to ^ when r^»£ and as shown in 

figure (2-0) below. This means that when M is not too small, 

b has a greater weighting in the determination of the value 

of r . so that an increase in the horizontal distance is 

o 

more effective in reducing the strength than the increase 
in vertical distance from the plate to the source. 
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(5) When the original sound source has a directivity 
(such as a quadrupole) , it is advisible to locate the slot 
in the direction of minimum radiation (see figure (2-P) below. 





\ 


\ 




\ 














3 


Figure (2-P); Direction of minimum radiation. 
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Chapter III; 

GENERAL METHOD FOR HIGH FREQUENCY WAVES; 

III-1 Background; 

The development given below indicates a general method 
that will be used in solving problems of acoustic wave 
propagation in the high frequency end of the wave spectrum. 
The method is developed out of a combination of the ray 
acoustics, Kirchhoff's integral and the stationary phase 

method and applied to the problem of a semi— inf inite solid 
plate. 

The basic principle can be outlined as follows. Consider 

a sinusoidal monopole source located above the semi-infinite 

plate. In the limit of very large wave number |</ the 

approximate solution can be constructed by means of ray 

acoustics. The solutions are defined in three regions with 

sharp discountinuities . The solution is reasonably accurate 

1 

and the order neglected is ^ . This solution is not 
accurate at the boundaries separating the three regions 
where a transition layer with contributions to the solution 
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of 0 ^ - are present. The proposed method uses the 

Kirchhoff s integral with a surface lying outside the 
transition layers. The solution obtained by ray acoustics 


defines the first approximate solution for (}) and 


an 


that are contained in the Kirchhoff 's integral. It is then 
possit)le to determine the function 4 > at any point p including 
that within the transition layer. This integral generates 
the second order correction which is of order 0(^). ..e 

method of stationary phase is used to determine the second 
order approximation directly by applying the stationary 


phase method to the integral, 


In a sense, the method is to treat the Kirchhoff 's 


integral as an integral equation which is solved by an 
iterative procedure using the ray acoustics solution as 


a first approximation. 


In order to demonstrate the development and application 


of the method, a three-dimensional problem is solved as in 
Section III-2 which follows. From the solution derived for 


the three-dimensional problem, the solution for a two 
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dimensional case is deduced and compared to already known 
solutions. Such a comparison shows a good agreement for 
the present solution with known solutions abtained by more 
elaborate mathematical procedures. 

An Illustrative Problem of Diffraction of Sound from 
a Three-Dimensional Source . 

A three-dimensional source S is situated above a 
semi-infinite plate as shown in figure (3-A) below. 

Problem Statement; The acoustic field due to a sinusoidal 
monopole source located above a semi-infinite plate in a 
stationary medium is to be determined. 
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Assumptions: (i) A simple (three-dimensional) source whose 

wavelength is short is used. 

(ii) A semi-infinite plane extends along the 
-oo<x^o region. 


(iii) The surface of the plane is rigid so 
that the normal velocity vanishes there. 

(iv) A stationary medium surrounding the source 
and the semi-infinite plane. 

Method of solution; With this arrangement, the first order 
solution can be obtained by use of the ray acoustics. Such a 
first approximation would provide solutions for regions I, 

II and III indicated in figure (3-A) as follows; 


Region I; 



^ikR gikR 

47TR 47VR 


Region II; 



ikR 

e 

4TTR 


1 


(3-1) 


Region III; J 

where R and R are the distances from the source and from the 
image to the observer at P respectively. 


This approximation is reasonably accurate as a first 
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solution approximation of the problem, but will not be 
sufficiently accurate in the thin shaded zones marking 
transitions from the source reflection reception to the 
source affected zone only; and that from the source affected 
zone to the shadow zone. 

The method now formulated will avoid these difficulties 
and hence provide a solution that can be applied at any point 
including the shaded (transition) zone. To do this, use is 
made of the Kirchhoff's integral as outlined below. Let G 
and <|> be defined inside a closed surface S as in figure 
(3-B) below. 



Figure (3-B): Boundaries for the Kirchhoff's integral. 
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The following equations can be written thus: 

= <5(y-‘X) (3-2) 

= o 

From equations (3-2) and (3-3) one obtains, by the standard 
procedure: 

4>v^<t - a v^4> = 4> <5 Cy-x) ,3-4, 


Hence 



r 

r 

c 

\ 

J j 



(J) V o 




(3-5) 


If the source considered is a point source, and the 
acoustic field is three-dimensional, then 



t=- 


Ae 


47T r; 


Q 


S 

4irr2 


(3-6) 


where ^ ^2^ 


X — X I as shown in the 
P ' 


figure (3-C) next page. 
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A. ik r 

4tL 



-7k K7 

e \7r; 


When k is large, and not too small, the second term is 
0(-^ as compared with the first term, and can be neglected. 
This gives: 

V c[) ^ Ajk vr: 

4Tnri 

vr = 

Ix-Xs| Ix-Xsl 


where e^_ is a unit vector from X to X ; so that: 
XS S 


CoS 6, — ^ 


I XS 


A. 

• n 


_ vr • n 


and thus: 


n . Vc|) cr - 


Aik o Cos 6i 

4*rrr, 


A similar evaluation can be done for the function Q; but 


in this case the normal vector n is negative and this gives; 


n. VQ =: 


-ikr, 


4irr2 


CoS e. 
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The Kirchhoff's integral equation (3-5) now becomes, on 


substitution of the above: 


rr 


=# -0r 

s 


^CoS0, + Cose^Js (3-7) 


with the assumptions of short wave length and hence k is 
large, use can be made of the stationary phase approximation 
in evaluating the integral appearing in equation (3-7) as 


shwon below. 


III-3 Solution by Stationary Phase Methos: 

The stationary phase method was developed by Lord 
Kelvin (21) toward the end of the last century and is based 
on the idea that the value of an integral containing a fast 
oscillating integral is essentially given by integrating 
over that part of the path where the phase change of the 
integrand is the slowest. The contribution from other parts 
of the integration parts will become very small due to 
strong cancellations. This method will be applied to the 


Kirchhoff's integral: 
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The surface of integration is shown in figure (3-B) . It 


consists of two surfaces bounding the left upper quadrant 


(surface S^)» and a spherical surface (surface 8j») thus 


completing the boundary. 


In equation (3-8) one may use the ray acoustic solution 

r dcf) 

for o, and r — . Those solutions are accurate to the order 

^ on 

of C ' except in the transition zone. The boundary shown 

in figure (3-B) is chosen to be outside the transition zone. 

/ 

Since tp and are known the itegral could, in 

principle, be integrated directly. Based on the assumptions 


stated previously, however, the mathematical procedure can 


be reduced immensely by the present method, 


The expressions for tj> and Q are, respectively; 

O 


•ikr. 

4TTKi 




ikr, 

e 

4'rrr2 


(3-9) 


For large k, the integral shown in equation (3-8) is 


approximately given by: 
S, Sr 


-ini 

dn J 


ik Cn+n) 

^ ds 


(3-10) 
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It will first be shown in the following discussion that 
for any given wave number k, the contribution to the 
integral from the spherical surface S|^ tends to zero as the 
sphere radius R tends to infinity. This result is obvious due 
to the fact that the distance from a point of integration on 
a spherical surface to the source S and the observation point 
P approach the same value as R becomes very large; so that 
the difference of the two terms ^ O as R— ^00 

^an dn ' 


Hence 




-8n 


dn J 


^ dS = o(l).0(|i).0(i|^)=oa): (3-11) 


n Ta 


o as R 


00 


From the ray acoustics solutions, it is observed that the 
integral that extends over the under surface of the plate is 
zero because this boundary is in the shadow region where the 
value of the velocity potential is identically zero. 

Equation (3-8) then reduces to the form: 



dS (3-12) 
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Figure (3-D): Half plane 


Let (X , Y , Z ) denote the position of the source 

S 3 S 


and (X , Y , Z ) the position of the observer and 
P P P 

j 77 j point of integration. The expression for 

r^ and t 2 appearing in equation (3-10) are defined by» 




cyp-'Vf + 


(3-13) U,b) 
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The method of stationary phase can now be applied if it 
is assumed k»l. For convenience, equation (3-10) is 
expressed in the form: 


® 00 

^(P) = ik jclT^j^ F (77,^) 


-00 





(3-14) 


where 




_J, 

/67T^ 



1 


(3-15) (a) 


and 

The function f = r^ + r 2 has a simple physical meaning. It 
represents the total distance an acoustic ray travels from 
the source point S to the observer point P through the point 
the integration path. As the integration point 
moves on the (TJ — plane, the phase angle k(r^ + r 2 ) 

of the exponential function in equation (3-10) varies; 
resulting in strong cancellations of the wave for large 
values of k. 


The attention will now be directed to the main 
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contribution of the integral in the neighbourhood of the 
stationary phase point (X) defined by: 


(tr;+tj) = o 


(r;+r^) = o 


(3-16) 


Obviously, ( 5 point that lies on the straight 

line joining the points S and P, as shown in figure (3-E) 
below. 



Figure (3-E) : Line joining the points S and P i 

passes through the stationary point. 


To apply the stationary point method, the function 
H-n,n is expanded into a double Taylor series around 
the point 
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where — — ^ t ) 

dV SC 


From equations (3-13) it is found that: 

i.(P; + r,')_ Itzk-hjzIL 
dTf' n r;. 

- !e=£ 


(3-18) (a, b) 




and 


dTj- 


(r+K'jrrfl+t') — 

(71-yS^ (yp-77f 1 

V ' 2; V17 Yi_) 

>i" r/ J 


STjdf^ 

.z 


(^;+>'2)=- 


(C-2.)(V-Ys) , fep-^Vy.-7?) 




Y-5 
'2 


(3-19) 
V (a,b,c) 


+rz) = (^,^ +y-[^£^ + _ 


At the stationary phase point Tj ^ Tj . ■• t- L , so that 

bV <7 Q 

and 


A 


4t('7 + »-,) 


)- Zp-t f 

17 1^ ) 


(3-20) (a,b) 


IT 

It follows from equations (3-13 a,b) that at ( 7^ • t) 

^ O 
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_ iV-Ysf (L-zJ 4 (VrvS Cep-. 


Xp 


(3-21) 


thus. 


xp-e Xp-4 


(3-22) 


That is, the three points ( ' ^s ^L' v^'O 

and ( XpjYpjZp") all lie on the same straight line. Using 

equations (3-13 a, b) , the expressions for ^(^1 “i* ^ 

( n + 7 — — (\~. + )^ ') may be written as: 

drjd^ ^ 


As 


|TI2/ = 


(n e) 


L. i_(2p-2s) 




inx 


(3-23) (a,b,c) 


where r - AKp C ^P“)s} subscript O 


denotes evaluation at 7^ zz and ^ ^ 



To carry out the integration separately with respect to 


and ^one can write equation (3-17) in the form: 




=: kr + + P' 


(3-24) (b) 


where the new variables 


(X^and p are defined by: 



-.2 







(3-35) (a,b) 


then 

} d^=J£ Jc-^" dn 

By substuting the expressions for A,B,C, one finds that: 



(Xp-Xs^) 

( r-zl 



and 
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The integral of <f> then becomes: 


6=jk [fx 

^ /67r^;Jr>z 


+2^£Ii£ 

Y] 


® ctna^ 


= jL- [fx. Xp.zXs J . 

/67T2 >Ai 7 r ^ 


-. 7*kr 
T_e 

J 


CX) 


i«\l IP 
e del 


dp 


'-C0 




(5-zG) 


where the end point value ^ will be defined later. Note that 
in obtaining equation (3-26) one has to use the relation that 


at the point U' 

^-Xs _ Xp-^ ^ Xp-Xs 

r 


•7 


^2 


(3-27) 


The definate integral in equation (3-26) can be readily 


integrated/ yielding: 


r 


CO 


10( 


' do(„ = 


- (1 + i) 


(3-28) 


—00 
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For the integral involving p one may define a function F(^) 
by: 


00 


Flf) = 



e’f dp 


(3-29) 


Upon substitution of equations (3-28) and (3-29) then 
equation (3-26) becomes: 

ikr 

<P =-^ F (p) 

4TT 


(3-30) 


If the common forms of Fresnel integrals are used: 


CCh) = 


o 


Cos ) dt 






SCZ)= sin 


J 


O 


F( p) niay be written as: 


F'P> = T-‘#[c(#P) + S(7fp) 


(3-31) 


(3-32) 


The asymptotic behaviour of C(Z) and S (Z) are given by the 


following equations. 
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C (s') acj + 




Cos 


7T 2 




(3-33) (a) 


TZ 


^ 2 — > ^ 2 > o) 

■*" o (-I3) 


(3-33) (b) 


It can then be easily seen that: 


ip 

F(p) c^l -6- _ 

^ )/W p 


F(P) ^1- 


1 + 1 


. e. __ 

yfW P 


+ 0(ia)i (F^°) 


(3-34) (a,b) 


In deriving the second expression in equation (3-34) (a) and 

(b) above, use has been made of the relations C(-Z) — -<^(Z) 
and S(-Z) = -S(“2) . 

The Lower Limit ^ : 

One recals that the equation (3-26) is integrated using 
the new variables defined by: 




(3-35) (a,b) 
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The line C( = O is the locus of point on which 

This means that along such a line, the function f has 

stationary values with respect to the variable ^ . The 

range of the limit of integration for p is, therefore, from 
* 77 

P to p oo as shown in figure (3-P) below. 



Figure (3-F) ; The r^nge of the limit of integration 
for p • 


The point P lies on the trailing edge of the plate where 
= O. From equation (3-24) (b) it can be seen that 

p = kr; — kr 

where r^= r^ + rj at the point P*. 


(3-36) 



To derive the expression r^ > one may write r^ and r 2 
the form; 



( 3 - 37 ) 





( 3 - 37 ) 


where ^ and ^ are the projections of r^ and r 2 on the 
9 ^ plane figure (3-G) below. 





Figure (3-G); Acoustic lengths on the projection 
of the plane (ij ) § ) • 


At 


* 

P 1 




satisfied the equation; 




(3 


Using equations (3-37) (a) , (b) , this becomes; 


in 


(a) 


(b) 


-38) 
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which can be written in the form: 



At the point ^ ^ one has, therefore, the following 

relation: 






Equation (3-42) is presented graphically by figure (3-H) . 



Figure (3-H); Graphical sketch of acoustic length . 

in relation to the stationary point. 



(3-44) 
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Conclusions and Discussions; 

Based on the geometrical configuration shown in figure 
(3-1) the main results of the above analysis can be 
summarized. 



Figure (3-1); Wave reception at P due to three-dimensional 
source in the presence of a semi-infinite 
plate. 

The formulae that have been developed in the analysis take 
into account only the radiation from the source S. To 
obtain the complete solution, the contribution from the 
reflected waves from the upper surface, which can be calculated 
by the same procedure, must be included. 

The formulae for calculating the acoustic field by a 
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source at S are as follows; 

Let denoted the potential at the point P due to the 

source S in the absence of the plate. 

<(>_(P) = - e^^^ (3-45) 

4 TTr 


The velocity potential <t>(P) in the presence of the plate will 
be: 

<|>(P) = F ( "p ) (3-46) 

where F ( p ) is the diffraction factor. The function P ( p ) 
is defined by the integral: 


• /— 

p 


(3-47) (a) 


= + ' S p)] (3-47) (b) 


= 


o 

X 

z 


(3-47) (c) 


p — ► 00 

p — ° 

1 P— 

€ 

The following three cases for the values of p are considered- 


Case 1 ^= 0 : 


This corresponds to the case when the line drawn from 
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the observer point P to the source passes through the edge 
of the plate. In such a case, 

F( ^ - 1/2 and therefore 

^ * 1/2 4*q(P) (3-48) 

This means that half of the radiation emitted by the source 
is received by the observer. 

Case 2 ^ < O; 

In this case, the observer is on the upper side of the 
line jB = 0. The edge of the plate lies below the 
stationary phase point. The diffration factor F( ^ ) is 
greater than half, and increases to unity, as the observation 
point moves further upwards. In the limit, takes the 

value so that the sound emitted is not obstructed. 

Case 3 P>0; 

This occurs when the observer is at the lower side of 
the ^ =s O line. The value of p is positive and continues 
to increase as the observation point P moves downward. The 
value P( p ) decreases as p> increases, and in the limit <|» 
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tends to zero, so that the observer is in the shadow region. 

For any given value of ^ , the variation of the 
diffraction factor P( ^ ) is shown in figure (3-J) below. 

As the observer moves from in the transmission region to 
P in the shadow region, the sound varies from complete 
reception to zero. The transition of F( jB ) from zero to 
unity occurs in the range of approximately (ref. 1) 

-20 ^ ^ 20 



Figure (3-J): Variation of the diffraction factor F(^). 
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One may note that for large values of k, the value of 
(j>(P) depends on the integration over a small region around 
a stationary phase point. If the edge of the plate lies 
inside this region^ it will make no effect on the transmission 
of sound from the source to the observer. 

Applications and Comparisons with known Solutions: 


Cooke (9) has derived an approximate expression for the 
far field due to a three-dimensional source over a 
semi-infinite plate. He gives the solution for for the 
direct radiation from the source as: 



Inifccp) + 1 S(f) 

- 2 . L 


(3-49) 


The solution derived in the present study is given by the 
multiplication of the results obtained in equations (3-46) 
and (3-47) (b) , so that 



(3-50) 


Apart from the difference in the factor ( — ^ \ ^ the 

4T 

solution obtained in the present study agrees with that 
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obtained by Cooke. 

The two-dimensional problem shown in figure (3-K) below 
may be considered as a special case of the above 
three-dimensional problem^ and it's solution may be obtained 
by letting Z s = O and P—^co in equations (3-44) and 
(3-46) . In taking the limit# the equations 

4>(p)= F (p) . ^ 

remain unchanged. The two important quantities that need 
to be calculated whenyO-^co are r, the distance from the point 
P to the source, and r’ojthe acoustic length given by equation 
(3-43) . 

n t 



Figure (3-K): Two-dimensional problem sketch. 
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In terms of the notations used in figure (3-K) , one finds: 


and 


•' + (@-e) 

^ ^-/OCos(©-e) as g 


00 


Using these values, the expressions of and ^ are 


simplified to the forms: 


ik/l -ikpCoS(9-e)^ 


<k(P) = 


4TT/S 

^-ikp CoS C®— 


> 


(3-44) 


(p 

(after removing the non-essential constant multiplier by 


re-normalization) and: 


p k (r;-r)= A/>( ii-G)s(®-8)) 


Cos (®^) 


Note that the sign of ^ is either negative or positive 
depending on whether the plate tip (^ = p ) is below or 
above the stationary phase point ^ = 0, figure (3-L). shown 
on the next page. For the incident wave, ^ is below the 
point B — O, so that it*s value is negative and is given by: 


?=-Ap 
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For the reflected waves, ® must be changed to < — ® and the 
value of B is positive as given by: 


p =Jzi^ CoS (V) 


The total velocity potential due to both waves is then: 
~7/c/oCos(0-®) 





(3-45) 




This result is identical to that given by equation (2.86) on 

* 

page (73) of Noble . When 8 is sufficiently far from ®— '7T 
and TT — © / the asymptotic expression for (j>^ for large values 
of kp may be obtained using the asymptotic expression of F 
given by equations (3-34) (a,b) . With the aid of the 

trigonometric identities* J2, CoS s: 1 *t* CoS (6±(B)^ 

2 Cos^^±S)as Cbsa + Cos© 


* Noble B., "Methods based on the Wiener-Hopf technique" 


(1958); where f> is denoted by r and F(U ) is replaced by 


-7T/4 

■ p— F('^) with FC^) = r 


e / dp 
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it is found that 


-ikyO CoS (e-®^ -ik/o CoS _ 

({> e + e + 


where 


ct £= (^Jkp) 


Z Sm ^Sio^ —ikfi+ 


Cos 0 + Cos ® 


<f>^ is the asymptotic expansion of the diffracted wave and is 


the same expression expressed by Candel (5) except for the 


sign in the complex exponent. 
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Chapter IV 

TWO-STREAM-FLOW ACOUSTIC PROBLEM 
IV- 1 Nature of the Problem i 

In conventional and short-take off and landing aircraft, 
the high velocity exhaust gas from the engine is usually the 
dominating noise source. The so-called "mixing noise", which 
is produced by the turbulent eddies in a region near the 
nozzle exit, together with the turbine and other engine noise 
is transmitted through a jet stream into the ambient 
atmosphere, which, in most cases, is in relative motion with 
respect to the aircraft. Figure (4-A) below shows the 
situation when the engine is installed over the wing. 



Figure (4-A); Over the wing external blowing. 
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Since the shielding of the wing is more effective on the 
high frequency end of the noise spectrum, the effect of 
shielding is studied based on the model problem shown in 
figure (4-B) below. 




Figure (4-B); Two stream flow with plane wave 
and solid plate. 


The problem consists of a wing surface which separates 
two of different velocities. In the region above the wing, 
the flow represents the exhaust of jet. The source of the 
noise is embedded in this flow at some distance from the 
surface. The flow below the surface represents the forward 
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motion of the aircraft, and the interest here is in studying 
the acoustic field below the wing. 

The following assumptions will be used for the study of 
this problem. 

(i) The problem is treated as two dimensional. 

(ii) The wing surface is approximated by a semi-infinite 
thin plate lying along the negative X-axis. 

(iii) The dorminant frequency of the source is high 
enough so that the short wave length approximation 
may be used. 

(iv) To simplify the mathematical treatment, the 
incoming sound wave is treated as a plane wave as 
an approximation to the case when the source is 
very far from the edge of the plate. 

(v) The jet stream is attached to the wing surface and 
extends over the entire half plane. 

The method to be used in studying this problem is 
basically that discussed in chapter III. The first 


approximation of the problem is obtained by applying the 
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method of ray acoustics. The boundanry surface consists of 
not only solid surface but also the interface between the two 
streams where transmission and reflection take place. 

In solving this classical problem (transmission of sound 
through a shear layer) , no consideration will be attempted 
in the stability of the problem. The study here will be 
confined to the study of the sound field below the plate, so 
that the boundary chosen for the Kirchhoff's integral consists 
of the half circle bounded by the flat plate and the free 
layer . 

The solution for this problem is obtained for the case 
of a large wave number. Since the problem involves no 
characteristic length, this solution will be valid for any 
value of k at a large distance also. The effects of the free 
stream, the deflection and reflection at the interface are 
contained in the results obtained. However, no numerical 
calculations have been attempted. 

IV-2 Problem Formulation: 


The problem outlined in section IV-1 above is 
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represented in figure (4-B) . The field equations are 
appropriately given as: 


^ -± (A. + UAjf = o 

dx^ dy^ 62 ^ Qo *6x/ * 


(4-1) 


^ I c j 1 1 <J f(i _ ^ 
-5-biF+ ° 


• 7r^ ^ ~ a: m- ~ 

where and <f >2 are velocity potentials above and below the 
y = O plane respectively. The boundary conditions are: 

O on y — O; -oo < X < O 
and = >^ 2 ^ on y = 0; O < X^i^oo where 


(i) 

iii - _ 


ay ay 

(ii) 

P.J = ?2 and 


and P 2 are the 


displacements above and below the y = O, X>0 plane 
calculated from and <j >2 

(iii) The Sommerfeld radiation condition is to be satisfied. 


The above problem will be solved by applying the 
Kirchhoff's integral over the lower half plane. To do this 
one has to find the boundary value of <f >2 at the interface by 
the analysis given in section IV-3 . 
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IV-3 (1) The Dispersion Relation of a Plane wave in a 
Uniform Stream: 


The convected wave equation in a stream of velocity U 


is 


d <f> <fi ± ( d , 
ax^ ay^ o^\et 



o 


(4-3) 


, 1 (k,x+l<2V)-lc*^t 

If € 


then equation (4-3) becomes. 


with the definition of 





(4-4) 


The locus of relation (4-4) is an ellipse as shown in figure 
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It is to be noted that for any given value of , two 
solutionns of K 2 are possible. The one denoted by"K^ has a 
possitive K 2 and corresponds to a plane wave propagating 
upward away from the y » O plane; while the one denoted by 
has a negative value of K 2 and corresponds to a plane wave 
propagating downward into the half-plane. 

^2) Reflection and deflection at the interface of two uniform 
flows with different velocities and U 2 : 

If two flows are involved as in section IV-2, then there 
will be equations for the two conic sections (ellipses) 
determined in the manner outlined in the last section. The 
two ellipses which intersect one another are shown in figure 
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The equations are given below: 


( 


W X X kz _ - 


Y'> o 


(4-5) 




kl _ k! 




/ 


(4-6) 


At the interface y = O separating the two streams and U 2 
an incident wave will produce a transmitted wave as well as a 
reflected one. As will be shown in the next section, the 
component of these three waves must be equal. With a given 
incident wave propagating towards the interface, the 
reflected wave K^. in the upper half plane Y > O is determined 
from the ellipse M = which is the Mach number of the upper 
stream, and the transmitted wave in the lower half-plane 
Y<0 is determined from the other ellipse denoted by M = M 2 
which is the Mach number of the lower stream. Thus, in the 
case M^ > M 2 as shown in figure (4-D) , sound waves are 
deflected toward the shadow region under the plate due to the 


presence of the velocity discontinuity. 
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(3) The Transmission and Reflection Coefficients t 

In this case the solution of equations (4-1) and (4-2) 
in the absence of the plate (that is, if the interface y = O 
extends from X = -oo toX = +oo ) are considered. The results 
are used to determine the boundary condition at Y = O, X > O 
in the Kirchhoff's integral to solve the two-stream problem. 

In the region y > O, the solution of equation (4-1) 
consists of an incident and a reflected wave while in the 
region y <0 the solution of equation (4-2) consists of 
transmitted wave only. If one selects: 


where 




-t 

7 “r 

<4 


= A; e 


= A, e 


1 

i (k,X + kzV-cot) 


4> = 


Ai. e 


j (k.xi-k'y-^t) 


(4-7) 


(4-8) 


where , A^ are the wave amplitudes for the incidence, 

reflection and transmission. 

The boundary conditions on the interface ( y ss O ) 



require that: 


(i) There must be a pressure balance across the interface 


hence 

P,=Pz 



where 




and 




or 

( 



on y=0 (4-9) 

(ii) 

The particle displacement of 

the two streams should 


be the same at the interface. This means that 
where the expression of the interface is taken to be: 

g (x,y,t) = V ('•-I'') 

and the equation of g (x.y.t) is: 

^ ^ “U* V q = O (4-11) 

at ^ 

V-Jhen this equation is linearized and applied separately 
to the upper and the lower streams, the following kinematic 


conditions result. 
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On matching these conditions, one obtains: 




(4-12) 


On substituting expressions (4-7) into equations (4-9) 


and (4-12) and setting y = O in the result, each term in the 


resulting equations will contain the factor e 


ik,x-icot 


A necessary condition to obtain a consistent solution is, 
therefore, to choose the same for all three wave solutions 
^1 / » snd This is what has been used in writing 

equation (4-10) and in determining and using the ellipses 
in section IV-3 (2) . 

Since the wave number k^^^ of the incident wave is given 

as a boundary condition, the value k^, ]x .2 Bud'ko for all three 

expressions in equation (4-8) are determined. The two 

equalities resulting from the two conditions at the interface, 

namely, equations (4-9) and (4-12), serve to determine the 

A 

reflection coefficient ^ and the transmission 






coefficient T^- Let D. =s OJ + U-K. and D_ = co + hr, . 

• » I ^ Z t 

i 
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Equations (4-9) and (4-12) yield; 


D, ( A. + ) = D^A^ 


*-i“ *t > = 


( 4 - 13 ) 


then 


2 12 

^ 2^2 - ^ 2^1 

2 12 

K^Di + k ' d , 
2 2 2 1 


"PlP2^2 

2 1 2 
K2D2 + K2D, 


( 4 - 14 ) 


IV-4 Solution of the Two-Stream Problem: 


The interest here is in determining the wave transmission 


and refraction in the region under the wing (see figure 4-E) ; 


for this reason the Kirchhoff 's integral is applied using a 


contour enclosing the lower half plane, and composed of the 


semi-infinite plate, the interface and a semi-circle of 


infinite radius. In addition, use is made of the short wave 


length approximation (that is^ k »i ) to find the 

asymptotic solution of <|>2 valid to the order 0(-^) . 
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Two important approximations are then made possible. 

These are: 

(1) The method of stationary phase can be used to evaluate 
the Kirchhoff’s integral. As shown in chapter III, the point 
of stationary phase is located on the part of the control 
surface which lies between the. acoustic *source and the point 
of observation. Thus, in the present problem the point of 
stationary phase lies on the surface composed of the solid 
plate and the interface. Any contribution from integral over 
the semi-circular arc is of higher order and need not be 
considered. 

(2) As shown in figure (4-E) the region of the transition 

layers approach two narrow sectors as k becomes sufficiently 
large. Outside these two sectors, the solution given by the 
ray acoustics is accurate to the extent that the neglected 
error terms are of the order . The solution sought for 

the transmission transition layer is of lower order ( that 

A 

is 0( hence can be determined by the Kirchhoff *s 
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integral using the boundary value on the y = O axis given by 
the ray acoustics. 

With these in mind, one may reformulate the problem in 
a way that one can solve it by the general method outlined 
in chapter III. Since the interest in this problem is in the 
transmitted wave across the interface, the field equations 
will be developed for the lower flow regime. 

IV- 5 Field equations and solution procedure. 

g = O (4-15, 

where L- _ ^ 

^ “ Q 

The appropriate boundary conditions are: 

(i) At the interface, y— 0; X>0, the solution of (j>^ 
is the same as the one given by the ray acoustics, namely, the 


form expressed in equation (4-7) as discussed in section IV-3, 

, , 7 (k,x-f k^y-oj-fc) 

4>,= A^e 

with and likewise determined by the method outlined in 


that section. Thus: 




d<fl • A i ^ 

dy 


o 

X>o 


and 


(4-16) 



(ii) At the backside of the plate, (that is, in y<0. 


X<0), the sound is completely shielded according to the ray 


acoustics. So, 


6 = ^ = o y =0 , 

^ ay 


X<0 (4-17) 


Let 


(iii) The radiation condition is satisfied at y 

ikMX , 

-U/ 


= e 


This gives, on substitution into equation (4—3): 


a. 1 r , K - [1/ = O 

• , ,-2 ^ \/2 ~ / I . * 


djt 

dx^ ^ 1-/V]2 ay-2 


with the boundary conditions of the form: 


— oo , 


(4-18) 


(4-19) 


(i) V=/\^e ■, |r = 7^^k,e i y=o 


x>o 


JL = k. - 


Mzk 




where 

(ii) 

(iii) The radiation condition at y — -o® 


W ^ dJE. -ss o on VssO,X<o 

dy ’ 


(4-21) 


On introducing new variables: 

" 2 * — ^ — / ^2 

X = X; y = yy 1-M2 ; K =y^ , K2 = 


y 1 -m; 


The equation (4-19) reduces to the standard form: 
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with the boundary conditions: 

lilX . 

(i) vj^srA^e 5 4| 

t ay 


•at 

1 /\t ki e 


on 


v = o 

y > o 


(ii) 


on 


y = o , X o 


(4-23) 

(a,brc) 


(iii) The radiation condition is satisfied at y — — oo 
Equations (4-22) and (4-23) constitute a standard wave 
scattering problem which has been discussed thoroughly in 
chapter III. 


The solution of these equations is: 

T(ii)<+k;y) 

k) =/lte F(P) 


(4-24) 


where 

and 

with 


F(f>. if re- “ 


00 ,,Z 


du 




as shown in figure 


(4-F) next page. The final solution for 0 is obtained by a 
direct substitution: - ; 

i [kf -y k*+ 2k Mjk, - (| -m|) kf y ] 

= e p(p) 


(4-25) 


where j 0 > ® > (T are defined next page. 




Chapter V 


CONCLUSIONS: 

The methodology that has been designed and used in the 
previous chapter proves extreemely powerful in handling the 
problem of current interest in noise shielding by aircraft 
wing. Additional features and utility of the mathematical 
method formulated can be illustrated in the following 
aspects of noise shielding. 

V-1 Shielding by a Finite Plate: 

The solution method for a semi-infinite plate can be 
used to obtain the solution for a finite wing. The following 
discussion is based on the two-dimensional case. As shown in 
figure (5-A) next page, the potential (f> (P) at a point P 
below the semi-infinite plate is due to the direct radiation 
and given by: 

<^(P) = + (^^ F(p) (5-1) 

where F(p is the diffraction factor, and 



98 


tv 



Figure (5-A) ; Geometry defining the diffraction 

factor for the two-dimensional case. 
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When P is very far away, ^ is given approximately by; 

k [r. -n Cos e'3 = sin^ , 

or 


F—yikiF Sfn f' 


(5-3) 


in this formula, 0^ is the angle between the line S ^ and the 
line OS figure (5-B) . 



Figure (5-B) : Far field approximation of wave 
reception at P. 


When P is very far away, this angle becomes approximately 
equal to $ , which is the angle between the line OP and a 
line passing through the two points S and O. According to 
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ray acoustics, the line OP defines the boundary line between 
the propagation and the shadow regions. 

The potential at P, to the first order of accuracy, may 
be written as: 

<l!> F(p) 

with p = Sm I 

As shown in figure (5-C) next page^the transition from full 
transmission to complete blockage takes place inside a sector 
centred at the trailing edge of the plate. The angular width 
of the sector decreases as the wave number K increases. 

Figure (5-D) next page shows the acoustic field due to the 
reflected wave. Again for large K, the transition takes place 
within a narrow sector. 

The implication of the above observation is that for the 
case of large K, the solution of ^(P) for a finite plate can 
be obtained by combining two solutions: one solution 
corresponds to waves diffracted around the trailing edge and 
the other corresponds to waves diffracted around the leading 


edge. 
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It can be shown that, for a long plate, the resulting solution 
will not significantly violate the boundary conditions imposed 
on the problem (figure 5-E) below. Using this approach the 
diffracted field around a plate of finite length has been 
calculated numerically. Figures (5-F, 5-G and 5-H) next three 
pages show the results of such calculations. 



Figure (5-E) : Configuration of a finite plate and a 
simple source above it. 


1. Incident and Reflected wave. 

2. Transition between 1 and 3. 

3. Incident wave only. 

4. Transition between 3 and 5 

5. Complete shadow. 
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V-2 Eff ect due to Gaps, Angular Deflection and other 

Geometrical Variation of the Wing Surface ; 

The audible range of the aircraft sound lies essentially 
between 20 to 20,000 Hz, corresponding to wave length of 
about 16 metres to 2 centimetres. If the dimension of any 
geometrical variation for the assumed plate and a gap on the 
plate surface is less than the dominant wave length of the 
sound, the fluid pressure and velocity variation over that 
part of the wing surface is nearly simultaneous. Therefore, 
the solution to the region of the acoustic field may be 
approximated by the solution of the steady compressible flow 
eguation. A local solution obtained by this equation is 
needed in order to smooth out any descontinuities in velocity 
pressure or mass flow due to ignoring such geometrical 
variation in the first approximation. In the gap problem 
discussed in chapter II, for example, the tangent velocity 
over the upper and the lower surface given by the wave 
equation and neglecting the presence of the gap will be 
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discontinuous and a local solution near the gap is necessary 
to remove this discontinuity. In the far field, this local 
steady solution around the gap acts like localized acoustic 
poles. The strength and phase of these induced poles depend 
on the nature of the original sound source. 

V-3 Effect of Source Distribution on Power Radiation ; 

The problems of acoustic interactions with barriers poses 
the more fundamental question of the nature the acoustic 
power is distributed in the region of interest. The clear 
understanding of this will facilitate the approach to the noise 
shielding. 

It is known that the interaction of acoustic waves by 
any barrier has the effect of inducing distributed sources 
on the surface of the barrier and hence the problem then 
becomes one of understanding the nature of interactions from 
such sources and the effect these have on the power generated 
from the original source. The study on such problems is 
nearly a century old beginning with the study of such problems 
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in the field of radio waves. Recently, however, the effect of 
a boundary on the radiated power has been studied by Thiessen 
(41), Embleton (12) and Levine (27), In such studies, it has 
been indicated that the distribution of sources affects the 
impedance of the power transmission. In this study, the 
radiative power of a three-dimensional monopole source placed 
in a stationary medium in the presence of a semi-infinite 
plate will be briefly discussed. The mathematical expressions 
given below are in terms of the spherical coordinates whose 
centre is the point on trailing edge that lies in the same 
plane as the source point S (figure 5-1) . The line through 
the point O into the paper is chosen as the polar axis whose 
azimuth and polar angles are 0 and 1 / respectively. 
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The ranges of such angles are: 


04 e^iir j o ^ T^/z 

The potential ^(p,t) at a point P(r,^ jQ) is obtained 


as in chapter III as: 

i (kr; -iot) 


<P=-M. 


4 irr; 


F(P,) - ^ 


i (kr^-ait) 


47 Tri 


F(Pz) (5-*) 


where r^ and r2 are the distances from the point P to the 
points S and S^- F( Pj") and F( p^) are the diffraction 
factors with and given by the expressions (assuming 


r very large) : 

p, = -JzkV^SmV Co3 (£^») ] 


z= +^Xkro Sini^ Gs 




(5-5) 




The function F may be written in the form F ss/[q *. 

I 

Aj( Pl) OfCPi) be denoted by and ; and 

by A2 ^nd ^2.* Then the following expressions for 


large r may be easily obtained: 

4irr 


e (5-6) 

4irr 



no 


I <p' > 

Po^O 


_ P^OO' 


32TT^a„r 


|A^+Ai+A,AiCos[^kCt;-ti)+ 

( 5 - 7 ) 


where/ for large r: ” ”“2l^ Sfn 6© Sirjl^ Sin 9 


The total radiated power W is then given by the expression: 


W = 


Po 60 ^ 

321T'"a. 



v4o 0=0 


A>Ai+A,A,, Cos[k(r,-r.')+(ar''i^ 


Smuc/pclB 


(5-8) 


which may be evaluated numerically. For the purpose of the 
present discussion, it suffices to comment on the following 
simple example. Consider the source to be at a distance h 
above the edge of the plate ( figure- 5-J) . 


X 


J£ 


I 

hi 

I 

I 

1 

1 


X 


Figure (5-J) 


Set up for power approximation for! 
distributed sources. , f 
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Assuming the wave number is very large, then the ray acoustic 
solution gives a good approximate solution. 

The acoustic field then comprises of three regions with 
sharp boundaries. In region I, which is the region of direct 
radiation and reflection, sound is received from both the 
original source S and its image so that: 


i(krx-wt) 

4TrPj 4ini 

The total radiated power from two identical sources with a 
separation distance 2/1 in a free space without boundaries can 
be easily be found. 



" 4TrQo 

where ^ 


[1 + i 


The power radiated in region I is 1/4 of this value and 


is hence equal to: 


W = . r + JL 

^ IblYQ, 2. LX J 


(5-10) (a) 


In region II, which is the region of direct radiation with 


no reflection, the radiated power is equal to 1/2 of the 
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It is seen that, as the source approaches the plate (that is 
H — ► o or |L( — ► O ) , the radiated power approaches a value 
which is 50% greater than the value for the source without 
the plate, so that 


Wt 

Wo 


J*S as /< 


On the other hand, when h is much larger than the 


acoustic wave length, then 

Wt _ 
Wo 


1 as p 


OO 


W-, 


Between these two limits, the value of 4 ^ oscillates 

Wo 

through a series of maxima and minima about the value 1 

Wo " 

The crossover points with 4^ -1 are M^nTT where n = i 2 — 

Wo ' J J 

The corresponding values of h is ( A = wave length) 

with y — ) 

As h is increased, ^ changes from greater than unit 

Wo 

1 to less than 1 as h passes the series of values h « 

n' 




The acoustic impedance as seen by the source is, therefore, 
greater or smaller than that in the absence of the plate 
depending on where the position of the source is. This effect 
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on the radiating power is due to the reflection effect of the 
plate. 

V-4 The Effect of the Uniform Flow ; 

Consider a monopole in a uniform flow in the presence 
of a semi-infinite plate as shown in figure (5-L) . 



Figure (5-L): A monopole in a uniform flow in the 
presence of a semi-infinite plate. 


The governing diffrential equation and the boundary 
conditions for the potential ^ are given by 


(5-12) 


k = 


Oo 


where 
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The boundary conditions are: 

fy = O at y = Oj X<0 

(ii) <f> and its derivatives are continuous elsewhere, 

(iii) The Sommerfeld radiation condition at infinity is 
satisfied. 

This problem can be reduced to the standard form by 
letting 


where: 



•j kMx 


f (5(,7,s,k) 


(5-13) 


x = x 5 y = y j 2 = 5 k = 


The equation and the boundary conditions for ^ are then: 

^ 2 ^ 

!|-z+ ^2. ^ ^ (5-14) 

(i) at y = O; X<0 

(ii) ijf and it derivatives are condinuous elsewhere, 

(iii) Sommerfeld radiation condition at infinity is 
satisfied. 

Therefore, the diffraction problem for in the new variables 
is identical to that for a monopole source in a stationary 


medium. 
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If the potential due to the incident wave only is 
considered the solution of is 

t=N;FCI) (5-15) 

1 - 

where W = — H and F( ft ) is the diffraction factor, 

^ 4ttP r 

Assuming, for simplicity, that both the source and the 
observation point lie on the X-Y plane, then ^ is given by the 
relation (figure 5-A)- 

B=k(s3 + Op - sp) 

In term of the original variables, the expression is 


P 


2 — 

= Jl. 

l-M^ 


yxo^+(i-M^)Vo^*-iyx ^+ -j (k-xof + ' 


(5-16) 


It will first be shown that along the straight line SOP^ 
shown in figure (5-M) the value of ^ is equal to zero for all 
values of K and M. Y ^ 


Figure 




(5-M); Geometry to show the boundary*; 
according to ray acoustics. ^ 
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Obviously, if SOP^ is a straight line, its end point (X^,Y^) 
and (X,Y) satisfy the relations; 


x=-rx, ) v = -Pv. 

where P is a positive constant. By a direct substitution. 


it is found that: 

a/x^+0-m')VJ +7x^4- y(X-X,f + 

(5-17) 


=yx?+o-M^)^^' + ryx'.t _y(|_p^)[x? tO-M"; v.’-J 


This proves that ^ has a zero value everywhere on the line 
SOP'^ for all values of the wave number k and flow Mach number 
M. When, in particular, the wave number becomes very large 
(K — ►CO ) ^ then the transition layer collapses onto this 
line SOP^ . Therefore, the line SOP^ is the boundary line 
between the regions of transmission and the shadow according 
to ray acoustics for all Mach numbers as long as the wave 
number k is large. Latter in this section, the effect of the 
Mach number on the transition layer thickness will be shown. 
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From this, it is concluded that the line drawn from the 
source to the trailing edge of the plate completely defines 
the position of the ray acoustics boundaries irrespective of 
whether the medium is stationary or moving figure (5-N) , In 
the special case of a plane wave incident on the plate, the 
above conclusion agrees with the result of Candel (5) as • 
shown below. 


The incident field due to a monopole at (r^ j 0^) and 
Z = O is given by the expression: 



Figure (5-N) : A monopole source field in two-dimensional 
96onietry in the presence of a uniform flow.i 
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When the source is moved to infinity along the line OS 
figure (5-N), the incident wave on the plate due to the 
concentrated source becomes identical to that of a plane wave 
with its potential given by the asymptotic expression of 
equation (5-18) with *-00 as; 




(5-19) 


The potential due to a plane wave of incidence angle ® 
figure (5-0) is given by the expression in the equation 19 of 


Candel (5) as: 



Figure (5-0): A plane wave field in two-dimensional 

geometry in the presence of a uniform flow.' 
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Comparing equation (5-19) and (5-20) it is seen that 
the incidence angle ® is equal to T\— *0o only when M — O. In 
general, when M ^ O, the incidence angle ©of the plane wave 
given by equation (5-20) is a function of 0^ and the . 
relationship is obtained by comparing the exponents contained 


in equations (5-19) and (5-20) thus; 


1 CoS0g "I _ O 

1-mL /l - J 1- 


Cos ® 


M CoS® 


(5-21)'(a) 


and 


sin 0o 


Sm ® 


(5-21) (b) 


vA — 1— MCbs0 

The following relation is then obtained from equation (5-21) 


(b) as; 


tan 0o = 


Sm® 


M “ CoS ® 


(5-22) 


This relation is also found to satisfy equation (5-21) (a). 

According to the conclusion mentioned above, the angle 
Oq is also the angle of the ray acoustic boundaries which 
were determined and denoted as < 5 + and 5 - by Candel (5)., 


figure (5-N) . Thus; 


ian 6-^ = tan (ir 0o') — — 


Sin ® 

M — CoS ® 


(5-23) 
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This relationship between and ® derived by the present 
general method is identical to equation 35 obtained by Candel 
(5) for the case of the plane wave. 

The effect of the Mach number on the thickness of the 
transition layer is now examined. An angle d is introduced/ 
for this purpose, using the line O as the baseline 
(figure 5-P) . The angle is defined as positive on the side 
adjacent to the radiation zone, and is negative on th£?^ shadow 
region. 

In terms of o( , the asymptotic expression of ^ for 
large r is obtained by a lengthy but straightforward 
calculation . 


|y 
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The final result is given as: 


B =_kjl^r A_M2sm=‘e.'+ M-^Sin (6o-«<) Sfoe,-a>s<^ 
b/ /l-M^ Sln=‘(eo-‘^) 


(5-24) 


(i) 

^ — ► — oo and F{ p ) 

(ii) 

P= O, E( p ) «l/2 

(iii) 

P -> + 0O , F(p ) 


+ O (’^Yr) 

This result is derived by assuming S and P are in the same 

plane. The sign of ^ is chosen so that ^ <T O when^V>0 and 

p t)0 when o( <C O • It follows from this choice that: 

->•1 as o( — ►oo 

cf=0 (5-25) 

->0 as — ►-00 

From the mathematical table of F( |5i ) (Ref. 1), it is 

found that F{ ^ ) is practically equal to 1 when p :Ss-20, 

The corresponding value of c( / denoted by 0^+ defines the 

position of the "bright" edge of the transition layer. 

Similarly/ when jB + 20/ F (P) is practically equal to 

zero, and the corresponding o( / denoted by o(— r defines 

the position of the "dark" edge of the transition layer. 

-5 , , 

As one can see from the expression for ^ / Oc^ and Of. 
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are functions of kr^ j 0^ and M. In table I and figure (5-Q) 
next page/ the values of and are listed. These 

values are calculated by assuming kr^ = 2500/ for 

various values of the Mach numbers. 

From this example/ it is seen that the thickness of the 
transition layer at the trailing edge of the plate becomes 
narrower when the Mach number increases/ and the decrease in 
thickness is more pronounced on the side adjacent to the 
shadow region (under the plate) . The effect of the Mach 
number is/ therefore/ to increase the size of the shadow 
region and hence make the shielding more effective. Table II 


shows such calculations. 
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Table I 


Transition layer thickness as a function of the Mach number . 


Mach number (M) 

0 

0.2 

0.4 

0.6 

0.8 

0.95 

0^^ (degrees) 

32.9 

IQII 

32.2 

31 .2 

29.5 

27.7 

(degrees) 

-32.9 

-32.0 

-29.5 

-25.5 

-20.3 

-15.7 

Layer thickness 
(degrees) 

65.8 

64.7 

61 .7 

56.7 

49.8 

43.4 
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Table II 

Values of ^ when kr^ = 2500 and * ^*5'^ 


0 < 

M = 0 

M = .2 

II 

• 

M r .6 

M = .8 

M = .95 

-35 

2 l ’.263 

t 

21 .820 





-30 

18.301 

18.759 

20.318 

23.763 



-25 

15.304 

15.667 

16.891 

19.546 

25.706 


-20 

12.278 

12.551 

13.463 

15.402 

19.670 


-15 

9.229 

9.419 

10.049 

1 1.358 

14.096 


-10 

6.162 

6.279 

6.660 

7.437 

8.982 


-5 

3.084 

3.137 

3.308 

3.649 

4.297 


0 

0 

-1 . 456 " 

- 1 . 019 “ 

- 1 . 651 ' 

- 2 . 543 ' 

- 3 . 863 "* 

5 

- 3.084 

- 3.125 

- 3.258 

- 3.512 

- 3.959 

- 4.522 

10 

- 6.162 

- 6.235 

- 6.464 

- 6.894 

- 7.628 

- 8.514 

15 

- 9.229 

- 9.322 

- 9.614 

- 10.154 

- 11.053 

- 12.105 

20 

- 12.278 

- 12.382 

- 12.708 

- 13.304 

- 14.277 

- 15.388 

25 

- 15.304 

- 15.411 

- 15.747 

- 16.355 

- 17.335 

- 18.439 

30 

- 18.301 

- 18.405 

- 18.732 

- 19.321 

- 20.264 

- 21 .307 

35 

- 21.263 

- 21 .361 

- 21.666 

- 22.216 
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